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We report the observation of two resonance-like structures in the tt^XcI invariant mass distribution near 4.1 GeV in 
exclusive B — > Ktt+XcI decays. A detailed Dalitz-plot analysis demonstrates that these structures cannot be produced 
by reflections from any known and possibly unknown resonances in the Ktt channel. If these two peaks are produced 
by resonances in the tt^XcI channel, their minimal quark structure would have to be a ccud tetraquark arrangement, 
similar to that proposed for the Z+(4430) structure reported by Belle last year in the 7r+ip mass distribution produced 
in B — > Kir+ip decays. In addition, we report some new measurements on the properties of the X(3872) meson and 
the 1 Y states that are produced with initial state radiation. The analyses are based on a large data sample 
recorded at the T(45) resonances and nearby continuum with the Belle detector at the KEKB asymmetric-energy 
e+e - collider. 

1. Introduction 

The £?-factories were built to study CP violation in £?-meson decays. They also serve as a plentiful source of 
cc — like quark pairs produced via the color-suppressed, but CKM-favored b — > ccs decay mode and by s-channel 
e + e _ — » cc annihilation via initial state radiation. A number of particles are observed, the so-called XYZ mesons, 
that do not fit into any of the predicted but as yet unobserved charmonium states. As a result, they are considered 
possible candidates for new, exotic types of particles such as multiquark states or cc-gluon hybrids. The former 
include either mesonic molecules (cqcq ) or diquark-diantiquark (cqcq ) tetraquark mesons. We report on the recently 
observed Z family of states that have non-zero charge, which means they can neither be cc nor hybrid states. In 
addition new measurements related to the X(3872) are mentioned, and properties of the 1 __ Y family of states 
that are produced via the initial state radiation process, including those of the possible ss analog, the 1^(2175), are 
reported. 



2. Charged resonance-like Z x and Z 2 states 

In 2007 Belle reported the first charged charmonium-like state, denoted by Z + (4430), in the 7r + ^ mass distribution 
produced in B — > Ktti/j decays [1]. After vetoing the peak regions of the if * (890) and if|(1430) resonances in order to 
minimize their effects (the if* veto), the tt+i/j mass distribution (shown in Fig. [1]) exhibits a strong enhancement near 
4.43 GeV/c 2 . A fit with an 5- wave Breit Wigner function for signal peak and a phase-spacelike function for the non- 
peaking background gives a mass and width of M = 4433±4±2 MeV/c 2 and T = 4bt\l(stat)tH(sysi) MeV/c 2 and a 
signal significance of 6.5<r. The product branching fraction is measured to be B(B° — > K ~ Z + (4430)) x B(Z + (4430) — > 
tt+t//) = (4.1 ± 1.0 ± 1.4) x 10- 5 . 

This Z(4430) + observation motivated a search for other charged charmonium-like states in 5° — » K~ir + cc de- 
cays 

The left panel of Fig. [2] shows a Dalitz plot of M 2 (if + 7r~) (horizontal) versus M 2 (7r + Xci) (vertical) for B° — > 
if _ 7r + Xci events. Some clustering is evident in vertical bands at the mass of the if* (890) and if | (1430), which 
correspond to B — > K*x& decays. In addition, a clustering of events in a broad horizontal band is also apparent. 
To characterize this signal clustering in the Dalitz plot, Belle did a likelihood fit with a two- variable (M 2 (if ~7r + ) & 
M 2 (7r + Xci)) function that is a superposition of various interfering Breit- Wigner amplitudes that includes all known 
if ~7T + resonances plus 7r + Xci resonances. The binned fit to the Dalitz plot is performed using a very small bin size. 



34 th International Conference on High Energy Physics, Philadelphia, 2008 




3.8 4.05 4.3 4.55 4.8 

M(nV) (GeV) 



Figure 1: The M(tt + ^ ) distribution produced in B — > K^i/j decays. The K* veto was applied. The shaded histogram 
shows the combinatoric background level inferred from side-band data. 




Figure 2: a) The M 2 (K~tv + ) (horizontal) vs M 2 (ty + Xci) Dalitz plot distribution for the B — ► Kn + Xci candidate events, b) 
The M(n + Xci) distribution projected from the second vertical slice in the range of 1.0 GeV 2 /c 4 < M 2 (K~7r + ) < 1.75 GeV 2 /c 4 
in the Dalitz plot a), which is the most sensitive to Z\ and Zi. The dots with error bars represent data, the solid (dashed) 
histogram is the Dalitz plot fit model with all known K* and two (without any) 7r + Xci resonances. 

To display the fit results, the Dalitz plot is divided into the four vertical slices shown in Fig. [2] a) and three 
horizontal slices. The second vertical slice in Fig. [2] a) is the most sensitive to possible Z + signals; the first vertical 
slice has some Z + sensitivity. The right plot in Fig. [2] b) shows the projected 7r + Xci mass distribution from the 
second vertical Kir mass slice, the region with most discrimination; the histograms are the results of fits to the entire 
Dalitz plot that are described below. 

The dashed histogram shows the results of a fit to the 7r + Xci structure with a fit function that includes all known 
Kir resonances, but no 7r + Xci resonance. The fit CL is very poor, only 4 x 10 -12 . When we include a non-resonant 
KttXcI term and an additional K* resonance with floating mass and width into the fit function, the fit CL (from 
a toy MC study) is still only 2 x 10 -4 . If we add one spin=0 7rxd resonance, the fit gives a mass and width of 
M = 4150±?J MeV/c 2 and T = 352±|| MeV/c 2 with signal significance of 10.7 a. However, the fit CL is only 0.1%. 

If we instead use two nxd resonances, denoted by Z\ and Z2, the fit gives masses and widths of Mz 1 = 
4051 ± 14t 2 ?MeV/c 2 , T Zl = 82±l\(stat)±il(sysi) MeV/c 2 , M Z2 = 4248t^(stat)t^°(^)MeV/c 2 and T Z2 = 
1 77I39 (stat)lg]; 6 (syst) MeV/c 2 , and a fit CL of 40%. The corresponding product branching fractions are B(B° — > 
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Figure 3: a) The M(J/^tv + 7i~) distribution for the X(3872) region produced in B° — > J/ipn^TT K^tt decays, b) The Ktt 
mass distribution for B — > X(3872)i^7r candidate events. The dash-dotted curve shows the non-resonant .KV component, the 
dotted curve the K* (892)° component and the dashed curve the background from sideband. 

K~Z+) x B(Z? tt+Xci) = (3Atli(stat)tl' 7 6 (syst)) x 10" 5 and B(B° K~ Z+) x B(Z+ tt+Xci) = 
(4.01q J^tat)!^ 7 (syst)) x 10 -5 . These are comparable to those of other char monium- like states as well as the 
Z(4430) + . The two resonances structure is favored over the one resonance assumption at the 5.7 a level. 

We tried to fit with many different Ktt resonance options including the addition of new resonances with floating 
masses and widths. However, no model that only has dynamics in the Ktt channel can fit the Dalitz plot. 

3. News on theX(3872) 

The threshold mass of the D°D*° system is 3871.8 ± 0.4 MeV/c 2 , while the PDG 2008 world average mass of the 
X(3872) is 3871.9 ±0.5 MeV/c 2 . This similarity in mass is a frequently cited characteristic of the X(3872) and highly 
suggestive of a DD* molecule-like structure. This motivated searches for X(3872) — > DD* decays, and peaks with 
masses that are slightly higher than 3872 MeV/c 2 were reported both by Belle (3875.2 ± 0.7±?;g ± 0.8MeV/c 2 ) and 
BaBar (3875.1±g;£ ± 0.5MeV/c 2 ). Both goups also measured branching fractions that are 10 times that for the decay 
B -> i\~X(3872); X(3872) -> J/^tt + tt~ . This might be trouble for the DD* mesonic-molecule interpretation Q of 
the X(3872), for which the predicted BF ratio is more like 0.1. 

Here we report results from a three-times-larger data sample of 657 million BB pairs and uses both D* decays to 
£>°7 and D°7r° to see whether the threshold DD* mass peak is really at a higher mass than 3872 MeV/c 2 . We also 
fit the signal with a more sophisticated form that uses a Breit Wigner and the Flatte fitting function. The Flatte 
function is essential if the resonance peak mass is below the DD* mass threshold. Fits assuming that the peak is below 
or above the threshold give the same signal significance of 8.8<r and a fitted peak mass of the Breit Wigner function 
of 3872. 6^0 4 ± 0.4MeV/c 2 , which is consistent with the current world averaged value for the X(3872) mass from the 
J/^tttt channel. We also measure the product branching fraction for the DD* mode to be (0.73 ±0.17 ±0.13) x 10 -4 . 

In addition, we searched for the X(3872) in the J/^7r + 7r _ mass distribution for the B° — > J/i/jtt+tt~ K + tt~ 
decays [6|. Figure [3] a) shows a 5a significance X(3872) signal in the J/^7r + 7r _ mass distribution from B° — > 
J/^7r + 7r _ K + tt~ decays and Fig. [3]b) shows the Ktt mass distribution for the X(3872) mass region. Here the non- 
resonant Ktt component, shown by the dash-dotted curve, is much stronger than the resonant if * (892)° portion shown 
by the dotted curve. The measured product branching fraction is B(B° — > X(3872)(i\~ + 7r - ) non _ res ) x B(X(3872) — > 
J/^tt+tt-) = (8.1 ± 2.0t^) x 10- 6 at the 90% CL and B(B° X(3872)(if *(890)°) x £(X(3872) J/^tt+tt-) < 
3.4 x 10 -6 . This feature is in contrast to corresponding values for conventional charmonium (cc) states such as J/i/j, ijj 
and Xd, where the B — > (cc)K* (890)° decay is much larger than non-resonant B — > (cc)K + tt~ contributions. 

We also report the first statistically significant observation of B° — > X (3872) Kg. The ratio of branching fractions 
is measured to be £ + Z^x ( 3 8 72 ) x+ ) = ^-82 ^ ^-22 ^ 0-05, which is consistent with unity. The mass difference between 
the X(3872) produced in B+ and B° decays is measured to be AM = M XK + ~ M XK ° = 0.18 ± 0.89 ± 0.26 MeV/c 2 , 
which is cosisitent with 0. 
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Figure 4: The M(A+A~) distribution with p tag in the e + e~ — ► A+A~ cross section via initial state radition. The solid curve 
shows the fit result of the signal. The dashed curve shows that of the threshold function. The combinatorial background 
parametrization is shown by the dashed-dotted curve. 




Figure 5: a) The e + e — > 07r + 7r cross section distribution fit with two coherent BW functions, b) The e + e — > </>/o(980) 
cross section distribution fit with one resonance and a coherent nonresonant term . 



4. Y states produced via ISR 

The Y states are 1 __ states produced via initial state radiation (ISR). Both the F(4260) [7] seen in J^7r + 7r _ 
masses and the F(4325) [9[ seen in ip masses were discovered in e + e~ — > JisrY processes by BaBar. They 
are not the same state. 

Belle confirmed both states with the same decay modes as BaBar's also using ISR. However, Belle observed one 
more peak in the TrnJip channel (denoted by F(4040)) with a mass of M = 4008 ± 40t^ 4 MeV/c 2 and the width 
of T = 226 ± 44 ± 87 MeV/c 2 with a significance of 7.4<j. This state's peak mass is consistent with the ^(4040) 
charmonium state in mass, but its measured width is much larger Belle also observed an additional peak in the 
channel (designated as F (4660)) with mass M = 4664 ± 11 ± 5 MeV/c 2 & width T = 48 ±15 ±10 MeV/c 2 , and 
a significance of 5.8cr [lOj. The F(4660) also does not fit neatly to any available charmonium assignments. Among 
these results, one interesting feature is that the Tnrifj' peak masses are different from those seen in the ttttJ/i/j channel. 

In addition, Belle observed a new near-threshold enhancement in the e + e~ — > A+A~ cross section using initial- 
state radiation [ll|. Figure [H shows a clear peak near the A+A~ invariant mass threshold. Here only one A c is 
reconstructed and the other is inferred kinematics when only the decay antiproton is detected. The signal function is 
a sum of a relativistic £-wave Breit-Wigner function plus a threshold function to account for non-resonant production. 
The fit result gives a mass and width of M = 4634± 8 7 (stat)±l(syst) MeV/c 2 and T = 92±f 4 (stat)±l° 1 (syst) MeV/c 2 , 
respectively. These values are consistent with the Y (4660) mass and width. This state can be interpreted as a dibaryon 



threshold effect such as the one first reported by the BES experiment in the M(pp) from the decay J/ip — > jpp [12]. 
However, the possibility of its being the 5 s Si charmonium state is not excluded. 

An interesting question is whether or not there exist any XYZ counterparts containing ss or bb quarks. Here 
we mention a possible candidate in ss sector. The F(2175) is a 1 __ state first reported by BaBar in the fo(j) 
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mass produced in e + e~ — » Jisr/o^ ISR events; the mass and width are M = 2170 ± 10 ± 15 MeV/c 2 and T = 
58 ± 16 ± 20 MeV/c 2 [l3[. Its producton and decay properties are similar as that of F(4260) and this has stimulated 
speculation that it may be its ss counterpart, the 1^(2175). This state was confirmed by the BES experiment in 
J/ij) — > r](j) /o(980) decays with a mass and width that agree well with the BaBar measurements fl4^ . 

We also measure the </>7r + 7r _ and (f)fo cross section produced from the e + e~ — > 7/s , J R</>7r + 7r _ and e + e _ — > jisR^fo 
initial state radiation events for the center-of-mass energy ranging from 1.5 to 3.0 GeV [15[. Figure [5] shows corre- 
sponding cross section distributions. 

A fit with two coherent BW functions is applied to the </>7r + 7r - cross section distribution; a fit with one BW 
function interfering with a non-resonant continuum function is applied to the 0/o(98O) distribution. By averaging 
the fit results near the mass of the Y"(2175), we obtain a mass and width of 2133t^ 5 MeV/c 2 and 169tg2 5 MeV/c 2 , 
respectively. The large uncertainty includes both statistical and systematic errors, where the systematic error includes 
the differnce of measured masses from </>7r + 7r - and (j)fo sample as well. This mass agrees with both the BaBar and 
BES measurements, albeit with a somewhat larger width. The fit to the low-mass peak in Fig. [5] a) gives a mass and 
width of M = 1687 ± 21MeV/c 2 and r = 212 ± 29MeV/c 2 , respectively, indicating that this is the 0(1680) radially 
excited ss state. The possibility remains that the F(2175) could be just a higher radially excited ss state. 



5. Summary 

We observe two more charmonium-like states, Z^ and Z^, with non-zero charge where both decay into Xcitt + . 
Their quark content is ccud. The X(3872) decaying into 7r + 7r - J/ip has been well established, while its underlying 
nature has not been conclusively identified. The favored interpretations are either a mesonic molecule or a diquark- 
diantiquark tetraquark meson. Some additional properties of the X(3872) are measured. The mass of the X(3872) 
decaying into the DD° final state is consistent with that from 7r + 7r - J/ip within errors. Another feature is that the 
non-resonant K+tt~ component in B° — > X(3872)if + 7r - ; X(3872) — > J/^7r + 7r _ decay dominates over the resonant 
if* (892)° component. This is not similar to the behaviour usually seen for conventional charmonium states. Among 
several 1 Y states produced via initial state radiation, we observe the F(2175) as a possible ss analogue to the 
Y"(4260) which, in turn, is a prime candidate for a cc-gluon hybrid. The F(4660) is also seen via its A+A~ decay 
mode. 
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